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Abstract 

Solitons are known to move ballistically through a medium without changement 
of their shape. In practice, the shape of moving inhomogeneous states changes, and 
a long lasting tail appears behind the soliton moving in a periodic medium. Such 
a behavior can be described within the model of exciton-phonon coherent states as 
a dynamic effect in soliton transport. We argue that the coupling between bosonic 
excitations of the medium, such as excitons, and elastic modes of it, such as phonons, 
can be responsible for these effects. We derive nonlinear dynamic equations for 
the excited medium in the long wavelength approximation (a generalized Zakharov 
system) and apply a kind of ballistic ansatz for the coherent state of bosons and 
displacement field. Like solitons, the quasistationary solution we obtain on this way 
can move through the medium ballistically. Unlike solitons and kinks, there are 
inhomogeneous corrections to the ballistic velocity and the coherent phase of the 
condensate that control the changement of the packet shape. In the limit of T — > 0, 
laws of conservation prescribe formation of the cloud of collective excitations around 
the quasistationary Bose-core of the packet, and the density of a growing boson- 
phonon tail behind the moving coherent part can be estimated. The total packet can 
be associated with an exciton-phonon comet with the quasistable coherent Bose-core 
and incoherent tail moving in the medium. 
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1 Introduction 



Some localized solutions of a special (and, actually, quite a narrow) class of nonlinear evo- 
lution equations are known to be able to move conserving their shape and interacting with 
each other like particles. The typical examples are the solitons of Nonlinear Schrodinger 
equation (NLS) and Korteveg-de-Vriese equation (KdV), the kinks of sine-Gordon equa- 
tion (SG) to name a few, [T[,@. It is worthy of mentioning that these particular solutions 
conserve all the intergals of motions prescribed by dynamic equations. However, as the 
evolution equations mentioned above are the result of approximation of more complicated 
dynamic equations, the solitonic properties of their solutions are the approximation as 
well. The question up to which extent they survive in 'reality' leads to interesting physical 
problems that urge people to go beyond the theory of exactly solvable models @,@],@. 

Fortunately, the exact results obtained within refine mathematical models can be 
compared with the real physical experiments. In addition, numerical simulations within 
more realistic models are performed for many equations supporting solitonic solutions 
Q. As a result, one can realize under which conditions the concept of soliton is a good 
approximation to interpret physical data. For example, the NLS equation is used to 
describe the phenomenon of Bose-Einstein condensation of dilute trapped gases and 
sharp pulses of light in optical fibers || , and the KdV equation is used in hydrodynamics 
to model the moving packets of surface waves in channels |§. The SG equation models 
surface grow and reconstruction in material science ]T0| . 

In this article, we discuss the transport properties of excited states (generally of elec- 
tron transition nature) in a medium (a periodic structure, such as a crystal or semi- 
conductor structure, some periodic biological tissues, etc.). Under certain conditions, a 
coherent state, which involves both the electronic (or spin) excitations and the medium 
elastic excitations, can appear in such structures. This state turns out to be localized 
inside the medium and can travel ballistically through it. In Biological Physics, such a 
state is known as Davydov soliton whereas in Condensed Matter Physics, one can 
mention, for example, Spin-Peierls Systems [12|, Charge Density Waves ||13|| , and anoma- 
lous transport of excitonic packets in semiconducting crystals and heterostructures |TJ| . 
In particular, we are motivated by the ballistic transport of excitons in pure 3D crystals, 
such as CU2O |T5| . However, the ballistic transport of localized coherent structures and a 
set of nonlinear equations being used to describe it can be found in modeling of different 
physical phenomena 0, [T(J, |T7[. For example, modeling of the coherent excited states 
of electronic origin ends up with NLS equation, whereas the coherent states of the elas- 
tic medium are described by the wave equation with nonlinear anharmonic terms (NLW 

In plasma physics, this is the case of Zakharov system of equa- 
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equation) 
tions |1| 



|,[23], while the Davey-Stewartson system can be mentioned in the context 
of fluid mechanics 
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Note that the two equations comprising such a system, NLS and 
NLW, are quite different from the mathematical point of view, but they are coupled with 
each other P^],P5|,[26] and this is the origin of the wealth of physical effects they can 
describe. 

Not surprisingly for the moving soliton-like states in periodic media, a shape of 
experimentally registered signals is far from being look like a "true" soliton (e.g., oc 
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l/cosh(LQ 1 (x — vt) ) ) or kink. One can assume that nonzero temperature, scattering on 
impurities and just a noise factor can be responsible for this. In some cases, however, 
the character of changement hints on some regular dynamic reason for these effects. For 
example, an asymmetric form of the signal with the pronounced sharp front and a long 



lasting tail behind the single soliton are observed on experiment [ 15 1. [27], [28 . Note that 
the subsonic ballistic transport is under consideration, i.e., (v) < c s , where (v) is the 
speed of the packet and c s is the (longitudinal) sound speed in the periodic medium. The 
similar effects, such as delocalized solitons, were also observed in different physical set- 
tings. For example, in a kind of pump and probe experiments with sharp pulses of light 
in fibers, the solitonic shape of a probe pulse was significantly changed by Stimulated 
Raman Scattering | 29| ; see also p0| . In fluid mechanics, one can mention the problem of 



wave breaking and appearance of the so-called peakons and compactons in the soliton- 
like solutions of generalized KdV equations |JT|; see also [j3"2] . In essence, the non-steady 
solitons are found to be able to move with an acceleration and generate "tails" behind 
themselves. 

In contrast, the transition to diffusion regime is possible for the ballistic exciton- 
phonon droplets, and, in the case of amplification of the soliton-like state formed inside 
the exciton-phonon packet, its ballistic velocity is found to be almost unchanged during 
such a non-stationary process |15|]. These facts hint at different physics controlling the 
moving coherent state formed by excitons and phonons. 

In this article, we assume that the coupling between collective excitations of the pe- 
riodic medium can be responsible for such effects in Condensed Matter Physics, and a 
proper dynamic model can reproduce them even at T — > 0. To support this hypothesis, 
we consider a field model that admits a soliton-like solution and can be easily generalized 
without a loss of physical clarity. 

This is a two field model, in which the excitations of the medium are modeled by 
a nonideal Bose-gas coupled with long wavelength modes of the displacement field. For 
example, two (or even several) interacting boson fields can model different phenomena in 
Condensed Matter Physics [fL8|l . In the case of semiconductors, one can often disregard the 
influence of free fermions or fermionic complexes on the processes under consideration. It 
turns out that many processes involving photons, excitons, and phonons can be described 
by use of the language of interacting Bose-fields p3 |. Moreover, if there is a branch of 



optically inactive excitons in a crystal, one can even exclude the photons from simple 
models dealing with such excitons. In addition, the lifetime of a moving exciton with 
hk ~ m x c s can be large enough in semiconductor structures, so that the transport 
properties of a packet of moving excitons can be observed experimentally ||15||,[p4 . 



2 Hamiltonian of the model 

The Hamiltonian of the medium can be taken in the following form: 

H = H gas (*P, V> f ) + H ph (u, vr) + H int (u, ^>). (i) 

Here, ip, if}^ are the Bose-gas field operators; they stands in the non-relativistic Hamilto- 
nian if g as , whereas u is the displacement field operator and tt is the momentum density 
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operator conjugate to it,. This pair stands in a Phonon Hamiltonian H pn . We consider a 
nonideal gas of phonons in the long-wavelength approximation and take into account the 
acoustic branch only, e.g., the medium is a 3D crystal of the volume V = LS±. 
The Bose-gas Hamiltonian has the following form: 



H gaa = Jd*E a ^(x) + — V^(x) + V ± (^( x )) 2 (^(x)) 2 + ^ (^t( x )) 3 (^(x 



where > is the strength of two particle repulsive interaction, and V\ > is the 
strength of three particle one. As the two particle interaction (~ u (^) 2 tp 2 ) can be 
strongly renormalized because of interaction with other fields, we include the hard core 
interaction term (modeled by repulsion ~ v\ {ip^) 3 ip 3 in H gas ). We assume that the 'bare' 
characteristic energies of the particle-particle interactions satisfy the following inequality: 

< Vi/cl^ < (<c) ^o/^x — const Ry x , const ~ 10, 

see 



for discussion. Here, a x and Ry x are the exciton Bohr radius and characteristic 
Rydberg energy, respectively. We count the energy of a free particle from E g = E g —Ry x > 
0, so that Ek = Eg + hk. 2 /2m. (For semiconducting materials, E g is the fundamental 
gap.) 

As only the longitudinal phonons U; interact with the Bose-field ip, ip 3 * in our model, 
one can exclude the transversal phonons (Vu ( = 0) from H ph . Then, it can be reduced 
to the following simple form: 

^^^^ pC? ^ 

— — h ~y d j u s d j u s (yi) + -y- k 3 d j u s d j u s d j u s (yi) + k a (djU s ) rfx, 



where q is the longitudinal sound speed of the crystal, and the dimensionless parameters 
k 3 and K4 account for cubic and quartic nonlinearities. (In this article, we will not take 
into account a quartic term, so k 4 ~ 0.) 

We take the gas-phonon interaction in the form of Deformation Potential: 

H int = J<t <9,%(x) ^(x) + tf a,-%(x) (-^ f A ^(x) ) dx, (2) 

where ctq and $0 are the coupling constants. Note that this Hamiltonian is equivalent to 

H mt = J O"0 djUj(-x) ■0^-0 (x) + 1? 9jUj(x) V-0^ V-0(x) dx. (3) 

Here, we choose cr > and do not fix the sign of i?o- Developing a theory, we have a 
freedom with the sign of $0; its value, however, can be roughly estimated as |$o| — fi 2 /2m. 
(Recall that <7 ~ E g .) 

If the exciton Bohr radius a x is no more than several times larger than the lattice 
constant a;, the second term in (|2[) can be important too. (In fact, this is an intermediate 
case between Frenkel exciton and Wannier one with a x ^> a;.) Then, the following coupling 
terms appear in Heisenberg equations 

E g %[) -> E$ + cr djUj 0, (4) 
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- (h 2 /2m) Aip -> -(h 2 /2m) Az/> - tf djUj Aip. (5) 

In terms of the lattice analog of the medium Hamiltonian, we made the hopping matrix 
element t (iy in H gas ~ £y of the 'lattice' boson ipj to be dependent on the 'lattice' 

deformation field iL 



Meanwhile, the energy on a cite, £ ,i — £ o ( m -^gas ~ £o,i$!%) depends on the lattice 
displacements too, f.ex., 



sec 



361 for discussion. 



3 Dynamic equations 

Our aim is to investigate a special class of solutions of Hamiltonian ([I]), namely, we 
will search for the localized moving excitations (or packets). In the case of $o — 0, 
such a solution exists and the packet can move through the crystal with the constant 
velocity v saving its shape in a manner as the solitons and kinks do in nonlinear media 
|IT|,P9|. To simplify the equations of motion, we choose the quasi-lD approximation. 
This means ?/>(x, i) — > ip(x,t), u(x, t) — ► u x (x,t). In fact, we assume that the packet is 
inhomogeneous along the Ox axis only, and v \\ Ox. Thus, it occupies all the area of S±. 



(For a discussion on the validity of ID approximation, see |37].) Then, we can write the 
Heisenberg equations of motion as 



ihd t ip(x, t) 



r, 2 



(E g -—<%-$o 9 x u x (x, t) dl + u ^^p(x, t)+u 1 ^ j2 ^ 2 (x, t) }i>(x, t)+a djUj(x, t) xj>(x, t), 

(6) 

d 2 -c 2 d 2 x )u x (x,t)~ 

- cf2K 3 dlu x d x u x (x,t) = p~ 1 cr d x {$'4>{x,t)) + p~ 1 ^d x (d x $ d x $). (7) 

After some energy (and momentum) was pumped into the medium, an excited state of 
it appears near the boundary where the external energy (and momentum) was absorbed, 
see Fig. 1. Due to initial conditions, the excited state is a localized one, and it can be 
modeled as a droplet consisting of excitons (Bose-particles) and phonons. Moreover, it 
can acquire an average momentum directed along the Ox axis because of unidirectional 
phonon production during the thermalization stage. Thus, the exciton-phonon droplet 
starts to move p9[ . The important assumption, however, is the appearance of a moving 
coherent field, or a condensate, from the localized excited state of a medium. For example, 
if the medium is a 3D crystal or an array of channels embedded into a matrix, it has to 
be coupled with a thermostat at low temperature. Moreover, the medium has to posses 
an inner mechanism of fast thermalization of the electronic excited states, for example, 



by emission of phonons [38 . 
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Figure 1: A medium, in which the boson-phonon soliton can propagate, is presented on 
Figures 1 (a) and (b) in the form of the channel 'abed'. It has the dimensions |ab| = L, 
| be | ~ VS±, and L 3> L ch , where L ch is the characteristic width of the soliton. After 
some amount of energy has been pumped into the medium during a short time interval St 
and absorbed near the boundary, a localized excited state is formed near the face 'ad'. It 
is schematically shown on Figure 1 (a) as a mixture of excitons and phonons. If there is 
a mechanism of the momentum transfer to the excited state, the droplet begins to move 
toward the opposite face 'be' with the velocity (v), see Fig. 1 (b). Such conditions can 
favor the appearance of a coherent boson-phonon state (an analog of Davydov soliton) 
moving ballistically along the axis Ox at T < T c . In other words, a sort of Bose-condensate 
can appear because of the effective attraction among the bosons (excitons) at T < T c , see 
Fig. 1 (b). The coherent state, however, is only a core of the total moving packet. The 
profile of the excitonic Bose-core, n (x, t) ~ \^/ (x,t)\ 2 , is shown by the dashed line and 
the intensity of the elastic (phonon) part of the Bose-core, d x u 0>x (x,t), is represented by 
changements of the intensity of the background color on Fig. 1 (b). 
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It is known [K0j that nonlinear elastic lattices support several types of localized excita- 



tions; some of them turn out to have the so-called nonzero direct current, u x ~ u (x — vt). 
For example, a moving packet can consist of two parts, such as 



Ut. ~ U n (X 



It) + J2C{x- vt) e^-^ b k + h.c, 



and v is the group velocity of the packet. On the other hand, the Bose-gas has a property 
of multiple occupancy, so that the coherent mode can be introduced through the substitute 
ip(x,t) — > ty Q (x,t) and / c?x |\E' (x, t)| 2 ^> 1. Thus, in the limit of T — > and within the 
long wavelength approximation, one can explore how Eqs. (H), (|7|) support the localized 
coherent solutions. In this article, we use the language and technique of Bose-Einstein 
condensation [[|l|] that are suitable to take into account both the effect of many-particle 



coherency and the effect of many-particle nonlinearity. Note that the two parts of the 
coherent field, ^ Q (x,t) and u Q (x,t), are taken into account selfconsistently and we do not 
integrate out the phonons. 

Substituting the coherent fields (c-functions) into operator equations (|6]) and (|7]), we 
obtain the following system of dynamic equations ($o = (h 2 /2m) $ and $o is dimen- 
sionless) : 

ihd t ty (x,t) = -(h 2 /2m)$ d x u o (x,t)dlty o (x,t) + 

h 2 

+ (E g -—d 2 x + v \* o \ 2 (x,t) + i / i\^o\ 4 (x,t)^ (x,t) + a d x u o (x, t) ^ (x,t), (8) 

d 2 u Q (x, t) — c 2 d 2 u Q (x, t) — c 2 2k 3 d 2 u d x u Q (x, t) = 

= p-V ^(^ o (x,t)) +p- 1 $ d x (d x %d x y o (x,t)). (9) 

This system can be considered as a generalization of the Zakharov system which appeared 
in Plasma Physics The main difference is that, in the case of collective phenomena 



in Condensed Matter Physics, Eqs. (H) and @ define the main part of the moving 
packet. (It is called the coherent Bose-core in this article.) In fact, in the case of T ^ 
and, generally, in all the nonstationary cases, these two equations have to be coupled 



with another system of equations on the so-called out-of-condensate excitations |IT] . The 
dynamics of these excitation states can strongly influence the dynamics of the "parent" 
condensate. 

In this article, we rely on a kind of adiabatic hypothesis to simplify the solution of 
the problem. If the occupancy of the coherent mode, N Q (t) = J dx |\l/ (x, t)| 2 , and the 
occupancy of the out-of-condensate cloud 5N(t) = J dyi{S^Sip(x,t)), change in time 
slowly and N Q (t) ^> SN(t), one can try to find some quasi-stationary solution for the 
coherent core ^ {x, t) and balance its energy and momentum by a kind of "leakage" from 
the Bose-core into the incoherent out-of-condensate cloud and tail, N — > N — 5N(t). 
This approach can be called the propagation of an exciton-phonon "comet" in a periodic 
medium, and it is valid within a finite time interval which is estimated below. 

Recall that the standard ballistic ansatz with v = const, 

u Q (x, t) — u (x — t v ) , 



7 



\l/ (x, t) = exp(i(ip c + k x)) exp(—iQ(k ) t) ip (x — tv), 

where 

h£l(k ) = Eg + (h 2 k 2 /2m) + /x, hk = mv, 

does not work properly in the case of t?o 7^ 0. Moreover, with such a ballistic ansatz, we 
always have v — v s , where v s>x oc d x ip c (x, t) = const is the superfluid velocity of the Bose- 
core. As the anomalous ballistic transport of exciton-phonon packets in semiconductor 
and quasi-lD structures resembles, to some extent, the well-known effect of superfluidity 
42] , we apply the methods of Many- Particle Physics to Eqs. (H) and (H). We assume 



that the following substitutes for the ballistic velocity and the chemical potential, 

v -> v{x, t)»»(l± ({v) ^ 2 (x, t) + ... ), (10) 

li -> /x(x, t) « /x (1 ± C(v) ^(x, £) + ...), (11) 

are more appropriate to model the nonlinear dynamics of the ballistic boson-phonon 
packets in the case of $o 7^ 0- However, within the model with Hamiltonian (^), one can 
start, for example, from 

v — >• v(x,t) « vT\ (d x u Q (x,t) ) 
and derive expansions (|H]) and (|TT). 



To simplify the dynamic equations, we choose the following ansatz: 

u D (x,t) = u (x - tvJ 7 1 (d x u (x } t) )J, (12) 

# o 0M) = exp(^ c + z k Q T 2 {d x u ) x - iVt(k Q , d x u Q )tj ip (x - t vT x { d x u Q (x, t) )J, (13) 
where 

ftft(fc , 9 x m ) « E'g + (h 2 kl/2m)J 7 3 (d x u (x,t) ) +/i J 7 4 (d x u (x, t) ). 
For JFj, we use the following expansion: 

Tj « 1 + Cj $ d x u (x, t) + Cj ( $ d x u (x, t) ) 2 , j = 1, 4, 



and hk = mv. Note that, in some sense, ansatz ( |13|) is the well-known representation 

*o(z, *) = \Jn (x,t) exp(i(p c (x, t) ), 

with the inhomogeneous fields n Q (x, i), u Sja; (a;, i) oc d x cp c (x,t), and fi(x,t) oc d t <p c (x,t). 
In this article, we will not write out and solve a set of nonlinear hydrodynamic equations 
f4"Tf on these variables. Instead, we are interested in a microscopic approach aimed to 
clarify the role of the coupling parameters, $o and ctq. Obviously, for a packet moving 
nonuniformly, one might expect some deviation from the simple low 

y? c (x, t) = tp c + p s x/ft - (Eg + p 2 /2m + /i) t/h. 

If a good initial approximation for such a deviation can be found, one can finish with 
a kind of relatively simple NLS equation on the envelope functions of the Bose-core, 



S 



exp(iif c ) ip (x) and d x u (x). Such a mathematical structure do appears from Eqs. (§) 
and @ within the validity of the adiabatic hypothesis if one starts from quasistationary 
ansatz (pD and fll3|). 

Recall that the following relation between i/) (x) and d x u Q (x) is a good approximation 
for the selfconsistent boson-phonon coherent state at i?o = 0> 

d x u Q {x,t) = f{^ 2 Q {x,t)) ^ A 2 tfj 2 Q {x,t) + A^ Q {x,t). 

In fact, this approximation is valid in a more general case (although it is a part of a more 
complicated relation, such as d x u Q = f(ip , d x ip ) ). Then, for the coherent phase <p c (x, t), 
we can write out the following symmetric representation: 

<p c (x, t) w <f c + k {1 + ( 2 (v) if%{x,t) H } x - 

- (E g + hk 2 /2m{l + ( 3 (v)tlj 2 o (x,t) + ---} - fi {1 + <U(y)tf(x,t) + ■■•}) t/h. (14) 
The ballistic velocity of the condensate, x — > x — v(x,t)t, is modified as follows: 

x -> x - v{l + d(v) ip 2 (x, t) + ■ ■ ■ } t. (15) 

Here, we introduce the parameters (j(v) ~ CjC( v ) ^0 an d expect (\(v) ~ (2(1^) — 
Cs(v) — div). If C2 — 1, we get immediately c\ = 1 + c 2 — 2 and c 3 = 1 + 2c 2 ~ 3 
within the quasistationary approximation for the dynamic equations. The next field 
terms ( oc Cj i?q i/)o(x, t) ) in expansions (|14D , ( |15|) can be defined through c 2 and c 2 as 
well. However, they are small in a weakly nonlinear case and we will disregard them. The 
constant c 4 will be defined to simplify the quasistationary equation on ip 2 (x). 

To clarify the dynamic properties of the generalized ballistic ansatz and define the 
unknown Cj, we will calculate the energy and momentum of a moving soliton-like solution 
with the factors J r j(x,t). Meanwhile, the envelope functions, exp(iip c ) ip (x) and d x u (x), 
are found from the quasistationary equations one has to derive. These equations read 

+c 4i '&od x u (x,t) H }ip Q (x,t) ps -(h 2 /2m)'dod x u (x,t) d 2 x ip Q {x,t) + 

h 2 

+ ( - l^fl + u o t) + v\ ipo( x , t) )ipo(x, t) + cr d x u (x, t) ip (x, t), (16) 

v 2 dlu Q (x, t) (l + ci ^0 d x u Q H ^ - cf dlu Q (x, t) - 

- c 2 (2k 3 d 2 x u x u o (x, t) + 3k 4 d 2 x u x (d x u (x, t)) 2 ) ps 

~ p-Vo^(V 2 (x,t))+p-^o«9 x (^:9^ (x,t)). (17) 

Note that we did not include the terms depending on time explicitly into these equations. 
The consequent divergencies, however, can be cured by taking into account the incoherent 
parts of the packet. In this article, we propose the effect of slow dynamic redistribution 
of the occupation numbers between the condensate and out-of-condensate cloud, while 
both of these states can be easily calculated in the quasistationary approximation. Recall 
that dynamic mean field interaction between the condensate and noncondensed part of 
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the packet leads to the damping of both of them [Q. On the other hand, the time 
interval, during which the exciton-phonon packet with the coherent Bose-core moves in 
the periodic medium, is always finite, especially for the "clean" samples without lattice 
imperfections. Therefore, the adiabatic hypothesis can be used in our case. 

To find a soliton-like solution of Eqs. (^) and (|17D , we reduce Eq. (|T7| ) to the following 
form: 

-(c 2 - v 2 ) d 2 x u (x, t) - c 2 2 R 3 d 2 x u Q d x u Q (x, t) = 

= p-^d^l&t)). (18) 

Note that the interaction vertices k 3 and o"o are slightly renormalized in this equation, 
^3 — %(^o)) ^0 — &o($o)- F° r example, we can write R 3 as follows (c^ — 1, c\ ~ 2) 

^3 = -|^s| + 2 (^ 2 /cf) |^o| < 0, (19) 

and \Rs\ < \ks\ for $ < 0. Here, we choose k 3 < (k 4 > 0) and assume the inequality 

|«s| > \M 

to be valid. Then, the sign and the order of value of k 3 remains the same after the 
renormalization. For example, we take |k 3 | ~ 5 — 10, whereas |$o| — 1 — 2, and, for 
v < ci, we can estimate v 2 /cf ~ 0.5 — 0.9. 

One can represent the solution of Eq. (|18D in the following form: 



d x u (x,t) « -A 2 ip 2 (x,t) +A A ^(x,t), 

where 



^0 



/p(c 2 - v 2 ) = j(v) (a /Mc 2 ) af > 0, 



A 4 \R 3 \ A\ > 0, 



and p ~ M/af, 7(f) = cf/(cf — v 2 ) > 1. Immediately, we can rewrite the factor T for 
the ballistic velocity ($0 < 0) as follows: 

T x {d x u {x,t) ) ^ T x (\^ {x,t)\ 2 ) ~ 

« \ + \c,M\€M-lM€M- (20) 

For $0 > 0, we estimate 

J-i(cUo(M)) -> J"i(|^ (x,t)| 2 ) « 

« l-|Ci(v)|Vo(^*) + Ci(«)^^)- (21) 
We use the following representation for | \l/ (a:, t)\ 2 : 

|* (x, t) | 2 = ^ 2 ( x - f)t ) -> $ 2 / 2 (a;/L ), 

where f(x/L ) is the dimensionless 'shape' function, |/(x/L )| < 1, $ D is the amplitude 
of the coherent state, and x — v(x, t)t — > x. The characteristic width of the condensate is 
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estimated as (3—6) L . Although the terms oc ipt i n Eqs. Q2"D] ) and (pi|) could be important 
in strong nonlinear regimes (as well as other omitted contributions, e.g., oc (d x ip ) 2 ), we 
assume 0(^)^0 < 1 an d GC^O^o ^ Ci( u ) i n this article. Then, we can roughly 
estimate the dynamic factors Tj as 

^■(|^ (x,t)| 2 ) ^l±c J ((v)$ 2 f 2 (x,t), Cj ~l, (22) 

and the question is how strong the factor T can deviate from 1. If |$ | ~ 1 and 
7(d) ~ 5 — 10 (v < Ci), we can estimate the dimensionless factor CjQ(v)/af that enters 
Eq. (HD as follows 

Cj CM/ a ! - c 5 \Ml{v) (vo/Mc 2 ) ~ 1. 

Then, the most important parameter in Eq. (^) remains to be af$ 2 < 1. 

Now, we can qualitatively describe the dynamics of bosonic core of the packet, see 
Fig. 2. The top of the soliton |\I/ ( x — v (x, t)t )| 2 — > $ 2 / 2 (x/L ) moves with 

v -> v ± (Ci(u) fto P (x/L ) ~ 1, 

whereas the slopes of |\I / (x,t)| 2 move with 

t; -> v ± lO-^Ci^) $ 2 ) v, U opc (x/L ) ~ 10" 1 . 

The ballistic velocity of the tails of |\l/ (a;, £)| 2 remains to be almost unchanged, i.e. = v. 
One can introduce the important parameter 5v t0 p as 

5Vtop ~ ±ci |C(«)| 

and the dimensionless ratio 

^ « ± Cl |C( W )|$ 2 ^/L . (23) 

Then, the time scale At, during which our quasistationary ansatz can be used to describe 
the dynamics of the core of the total solitonic packets, can be roughly estimated from 
|5v top | At / L Q ~ 1. In fact, parameter ( p^) depends on time because the dynamic leakage 
from the Bose-core into the tail and coma of the "comet" leads to $ (£) and L Q (t). To 
estimate the values of \fi\, af $ 2 , and L , we have to solve the stationary equation on 
ip (x). 



4 Solution of stationary equation 

As a result of all the simplifications we made, we obtain this equation in the following 
form (we choose C4 = 1): 



( - f-{ 1 + ((v) ^ (x) - C» rt(x) } d 2 x - \v (v)\^(x) + vi{v) rt(x) )Mx)- (24) 



h 2 
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This equation resembles combined Ablowitz-Ladik and Nonlinear Schrodinger Equations 
taken in the continuous limit [36fl,|43|. Note that the boson interaction vertices are 
strongly renormalized because of interaction with the coherent phonon field d x u a . We 
assume the following conditions to be valid 

Vq > -> v^{v/c h a ) < 0, v x > ->• v\{vjc h a , R 3 ) > 0. (25) 

If the dimensionless parameter ((v ) $^ oc af $^ is small enough, one can reduce Eq. 
(|24] ) to Nonlinear Schrodinger Equation. We use the representation 

M x ) = $ f(P($o)x, v($oi), n = -\i4 < o, 



with unknown functions /3($ ) = L 1 and r)(& ). Then, the following simple equation 
appears as an approximation of Eq. (|24[) : 



\^\f{x)-\v {v)\e o f{x) + v l {v)^ o f{x) « !tf(§ )dlf{x). (26) 

2m 

This equation is the so-called "subcritical" NLS equation. In Condensed Matter Physics, 
it is used in the theory of superfluidity and Bose-Einstein condensation JHj]. It is also 
applied in nonlinear optics for light pulses in the medium with a cubic-quintic nonlinearity 



45| , and it is known as the Lienard equation in the theory of exactly solvable nonlinear 



equations [fig] . 

Here, we apply it to describe the coherent state of excitons and phonons with macro- 
scopic occupancy. The following approximation is used for the renormalized vortices: 
uq ~ and v\{y) ~ v~\{y ) + |z4)(f)| C(v). To estimate their strength, we use the following 
formulas: 

\Vo(v)\ = j(v) (ao/Mcf) ((7 af) - u > 0, (27) 
which as valid at 7(f) > j , or, equivalently, v Q < v < q (for discussion, see JT5|] , PT[ ) , 

p x {v) = ( 7 ( V ) |k 3 | ){l{v) (^) } 2 (a af) + v u (28) 

and 

^^«^^ + N$ 2 (C(^)^). (29) 

Note that in the case of $0 < we obtain the enhanced three particle vertex, and the 
nonlinear corrections are important to estimate its value, V\ — > V\ — > V\ > 0. 

The localized solution of Eq. (|26"D exists if the following inequalities are valid | 46| : 

Then, = |/z|(<3> ) and we have 



12 



As a rough estimate, we can use the inverse formula $ 2 (fi) ~ 0.5 ( \ji\/ ji*) $* 2 valid at 
|/i|//i*<l. 

The representation = {h 2 /2m) L 2 leads to an easy estimate of the characteristic 
width of the soliton. Indeed, we can introduce the length through ft* = (h 2 /2m) L~ 2 
and obtain the following representation 

/?($„) - /?(//) = ^fi*(\fi\/fi*), (3{ \n\ ) = L- 1 sj\fi\/fi*. (31) 



Then, L = and, always, L > L*. To estimate the important microscopic 

parameters fi* and L*, we write the renormalized vertices in the following form: 

\uo(v)\ = a (v) Ry x a x and v x {v) = ai(v) Ry x a x , (32) 

and estimate a — <5i — 10 _1 . (In theory, a can vary from ~ 10~ 2 to ~ 1 with 
changement of v within v Q < v < q fl9j, ||3~7ll .) Then, we have 

/i» = 1 ^ Ry x ~ 10~ 2 Ry x , Ll(v) ~ ^ a 2 ~ 10 a 2 . (33) 
lb «i «q 

As the following inequality is valid 

we can define the upper limit of Sv top /v ~ Ci(' u ) i n our m odel. We have 

CiOO ^) « l*o| Cl | 7 («) ^} (af $: 2 ) * 0.1 - 1, (34) 

where 

a? $:» ~( fib(t;)/ai(t;) ) (a?/a») ~ 0.1. (35) 

As a result, one can define the meaning of the weakly nonlinear case in our model. For 
example, if the parameter $ 2 /$* 2 < 0.1 — 0.5, the assumption 

C(v) ®t « ((v) $ 2 

we made to simplify the factors Fj( d x u Q (x, t) ) to J 7 j( ip 2 (x, t) ) is correct. 
We write the exact solution of Eq. (EH) in the following form: 



^ {x) = (36) 

yjl ~ M/V cosh 2 (/?(|/i|) x) + (1/2)(1 - yjl - \ii\h?) 

Note that 

^ 2 (x = o) = $ 2 ( |Af | ) = $: 2 (1 - v/l-H/zi*). 

If \fi\ <C /i*, we can use the following asymptotics (r)(\fi\/ fi*) — > 0): 

^o(x) ~ 2$ exp(- / 3($ ) |x| ) at |x| > 2(3{§ )-\ (37) 
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In addition, the coherent phonon part u Q (x,t) can be rewritten as follows 
d xUo (x) » - 7 («) (ao/Mcf) (a? $: 2 ) (^/$ : 2 ) / 2 (/?(|/i|) s) + 

+ 7(«) l%l {7(«) (V^) (a? *: 2 )} 2 (^: 2 ) 2 / 4 (/3(|//|)x). (38) 

We apply the 3D normalization condition on the macroscopic wave function ty (x,t) 
as follows 

J\^ \ 2 {x,t) dx = N (t) > 1. 
Substituting the quasistationary ansatz into the integrand, we have 

S h 2 {x-v(l + Ci(v)i/; 2 (x,t) + ...)t)dx = N (t) -> (39) 

-> 5 !ij)l{x')dx' = N Q . (40) 

Here, iV ^> 1 has a meaning of the number of particles inside the quasistationary conden- 
sate, (it is a macroscopic number), and S ~ ~ S±, where S± is the cross section area 
of the crystal. Note that, due to the symmetry of f(x/L ) against x — > —x, the value of 
A" is conserved, but first approximation. 

To understand how the macroscopic wave function is formed from the microscopic 
parameters of the theory (they are fi*, $*, L*, and the renormalized vertices v(v)), it is 
useful to introduce the parameter N* = S $* 2 as well. We estimate N* as follows: 

N*(v) ~ (l/fiaxtvj) (S/4) » 1. (41) 

Then, N*(v) ~ (S/a 2 .) is always the macroscopic number in 3D case (oci(v) ~ 10 _1 — 
10~ 2 ). In fact, Eq. (|40| ) leads to an algebraic equation that relates the macroscopic 
parameters with the microscopic ones. As a result, J\^\JJjl* can be defined as a function 
of exp(2 N Q /N*). To avoid cumbersome formulas, we use simple estimates: 

y J\ fJ ,\/ fJ ,* = f{N /N*)~N /N* and L /L, ~ (NJN*)' 1 - (42) 

They are valid up to A" G /A^* ~ 0.1 - 0.5. 

Thus, we can estimate the important dimensionless parameter that characterizes the 
deviation of the ballistic velocity v(x, t) from v = const and the chemical potential /i(x, t) 
from = const. We can write 

®l/K 2 = f( Wlf ) « (l/2)(H//0 ^ (l/2)(iV /iV*) 2 . (43) 

This means 

$ 2 /$* 2 ~ 10~ 3 -> 10~ 2 while NjN*(v) ~ 10" 2 -> 10 -1 , 

and 

$ 2 /$* 2 ~ 0.2 -> 0.4 while N /N*(v) ~ 0.6 - 0.7 -> 1. 
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As a result, the reasonable values for the macroscopic parameter ((v) $^ ~ 8v tov /v can 
be estimated as 



the interplay between the microscopic and macroscopic parameters of the theory: 



On Fig. 2, we show how the coherent part of the packet moves in the medium if one 
switch on the interaction $o 7^ 0. The initial state is taken to be symmetric (the exact 
solution for do ^ 0, "#0 = is used), and its evolution is obtained by use of quasistationary 
ansatz (|12j ), ( |i~3|) which conserves the amplitude and the characteristic width of the moving 
soliton. Although it is a rough estimate of the real dynamics of the Bose-core, it gives some 
understanding of how the microscopic exciton-phonon interaction controls the dynamics. 
In addition, it can be easily adjusted to a more realistic case of the moving core and 
out-of-condensate cloud. 



(the microscopic C(t>)$* 2 





(45) 
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5 Energy of the moving soliton 

Recall that the standard ballistic boson-phonon soliton, *ff (x,t) = exp(icp c (x,t) ) o (x — 
vt) and d x u Q (x — vi) with ip c (x, t) = <p + k x — u>o t, is a stationary state (#o = 0), that 

18 E. = f*VM = o»* and Po ^ = / TO , t )=co„ S , 

At i? 7^ 0, we obtained an analog of such a solution, starting from Eqs. QI2D, flTB| ) 
with f — > t>(x,t), /c — >■ fcoO 3 ^)) an d A* — *■ n(x,t). In this article, the quasistationary 
approximation is the most important one we used to simplify the dynamic equations for 
the core of the packet. Therefore, the exact calculation of E Q leads to d t E Q ^ for 
\l/ (x,t) = exp(iip c ( x, t, ip (x, t) ) ) ip ( x — v(x, t)t) with N Q = const and L = const. 
Within the validity of the adiabatic hypothesis (slow exchange between N Q (t) and 5N(t) 
while N Q (t) ^> 5N(t) ), one has, first, to calculate d t E Q of the quasistationary Bose-core 
and, second, to cure the obtained divergency by d t E tot = d t E Q + d t SE = (T — > 0). 

To calculate the energy of the moving condensate, we have to integrate the zero com- 
ponent of the energy-momentum tensor over the spatial coordinates. We have the 
following formula (written in the laboratory frame): 

ft 2 

T °(x, t) = E g m*J> + — V*: V* + ^ (* G ) 2 * 2 + ^ (^o) 3 *l+ 

2 2 fc2 

+ ^(^o) 2 + ^(^ ) 2 + ^«:3(^ ) 3 + a d x u o %V o + A^d a u V%V^ . 
2 2 3 2m 

After integration of T®, we conclude that there are no terms cx t l in E (t), 

E = E + (6E) (vt/L ) 2 + ■ ■ ■ « e N a + 6e N (Sv top t / L ) 2 . (46) 

In fact, the dependence on N Q is a nonlinear one, 

E « e Q (N /N*) N Q + 5e (N o /Nj N Q { 5v top (N ) t / L (N o ) } 2 . 

We found 5e > for both $ > (6v top < 0) and $o < (6v top > 0). Indeed, we can 
write out the expansion 

5e « 5e{, 0) + Se^ (C(v) $ 2 ) + • ■ ■ ~ 5e { ° ] > 0, 

and^eQ 1 ' 1 ~ Se^\ Recall the structure of the stationary part, E Q = e$ N D = E^+E^+E^. 
We have 

e « (h 2 kl/2m)( 1 + const $ 2 + ...) - const' |/i| + e ph , (47) 

where const, const' ~ 1 and e p h ~ e p ^ — e p ^. Here, the harmonic contribution into -Eph 
is e {2 l N > 0, and 

e «^Ai(^j 7W ^y (o 3 $ , 2){05|/i|//i . } 

16 
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Figure 2: To model transport properties of the boson-phonon soliton in a periodic 
medium, we started from the symmetric soliton as an initial condition at t — 0. Dynamics 
of the boson (exciton) part of the packet is presented on Figs. 2 (a) and 2 (c) in form of 
the moving l^o^ — v (x,t)t )| 2 without any "leakage". The phonon part (a moving kink 
of the displacement field of the medium) is depicted on Fig. 2 (b). (As a rough estimate, 
d x u Q (x — v(x, t)t) ~ — A 2 ipo(x — v(x, t)t). ) The interaction parameter (±(v) &l ~ Sv top /v 
controls the changement in the packet shape. It is taken to be +0.1 (Figs. 2 (a) and (b)) 
and —0.1 (only the Boson part of the packet is presented on Fig. 2 (c)). Then, the visible 
changements occur after the packet has traveled the distance of nL corresponding to 
Sv t0 p At/L ~ 0.5 — 1. For 5v top /v = ±0.1, we estimate n ~ 10. Note that the presented 
result is a crude estimate of the dynamics of the Bose-core of the total moving packet. 
To proceed, one has to "dress" such a core with the out-of-condensate cloud and tail. 
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We can compare formula ( |7 ) with our estimate of SeQ 



5e ~ (h 2 kl/2m) {cl/icf) consti + 3const 2 + const 3 > 0, (48) 

where all the constants are small (as a rough estimate, const,,- ~ 10 _1 ). Then, 5cq <C e$. 

The result (^) means that changement of the shape of the Bose-core, which was 
included into the generalized ballistic ansatz to satisfy the dynamic equations with $0 7^ 0, 
costs energy for both 5v top > and 6v top < during the observation time. Note that the 
width of the soliton, (3 — 6)L (N o ), Lq 2 cx \/i\(N ), and the number of Bose-particles 
forming the soliton, N Q , are conserved as it was presented on Fig. 2, but this is only the 
first approximation within the adiabatic hypothesis. It is possible, however, to balance the 
energy of the total boson-phonon packet. (Here, the image of an exciton-phonon "comet" 
helps: the Bose-core with iV ^> 1 is only a part of the moving delocalized packet.) 

First, we can calculate the value of { d]y E Q } 5N. We take it as 

{d N0 E o }SN ~ r 5N, 



where we have the following formula for eo, (eo 7^ eo in Eq.(47) ): 

eH ~ {h 2 kl/2m)( 1 + const new C(«) $0 + •••) ~ const' new + e ph , (49) 

— (2) (3) 

and e p h ~ 3e ph — 5e ph . Second, we assume that there is a "leakage" from the moving 
localized state, and such a leakage occurs dynamically, 

5N(t) ~ -SN (t) ~ - const A^ (Sv top t / L f < 0. (50) 

Then, the tail behind the condensate grows as 

(SN xMil ) ^ 6N Q (t) = const' N Q (6v top t / L ) 2 > 0. (51) 

To simplify the model, we disregard the creation of inside excitations, or the out-of- 
condensate cloud around the Bose-core, in this article. Then, const ~ const'. This means 
the following terms, 

- eo const iV (5v tov t / L ) 2 + 5e N (Sv top t / L ) 2 (52) 

and 

E tail {t) ~ (hu x ) 5N x , tail (t) + (hu ph ) 5N phMil {t), (53) 

have to be included into the (total) energy of the moving packet. The time interval 
At, during which our quasistationary solution for the core + tail can be used to de- 
scribe the transport of the total packet, can be roughly estimated from the condition 
\Sv top \At/L ~ 1, (see Eq. (|45D and Fig. 3). To estimate the excitation energies 
(hujn) and (huj p h) in Eq. fl53"D, we have to discuss the excitation spectrum of the moving 
condensate. 
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6 Excitations and the tail of soliton 



Within the quasistationary approximation, the following decomposition of the field oper- 
ators is used [f7H)|41||: 

$(x,t) = ® {x,t) + <fy(x,t), 

w(x, t) = u Q (x, t) + <5u(x, t), 

where the fields 5ip and 5u(x, t) describe the out-of-condensate particles. To consider 
the case in which the moving condensate emits excitations during the observation time, 
one has to introduce the fluctuation of the condensate, 5^ Q (x, t) and 5u Q (x, t), and the 
fluctuation of the quasistationary out-of-condensate cloud described by 5ip' and 5u'y If 
the number of the "lost" particles, 5N, is small during the observation time (but 5N(t) 
is continuously growing), one can consider the moving coherent packet as a quasistable 
one. Therefore, we can write the field decomposition as follows: 

■0(x, t) = exp(iip c (x, t, tp (x, £)) ) {^ ( x - v(x, t)t ) + 5tjj ( x - v(x, t)t, xj_, t ) } + 

+ 6V {x,t) + 64/(x,t), (54) 

Uj(x, t) = U Q (x — v(x, t)t )5ij + 5Uqj{x — v(x, t)t, Xj_, t ) + 

+ 5u (x, t)5ij + 5u'j(-x, t). (55) 

In this article, we do not take into account the fluctuational parts in this decomposi- 
tion and, strictly speaking, we can only estimate the the kinetic effects, such as the 
proposed d t N (t) ^ and d t 5n(t) = d t {8ip^ Sip)(t) ^ 0. However, the linear quasista- 
tionary equations on 5ip Q and Su j can be used to find the excitation spectrum of the 
outside excitations that is required for Eq. (|53f). Then, assuming the quasistability of 
the Bose-core during the observation time At, one can switch on the mechanism of oc- 
cupancy redistribution between the macroscopically occupied mode "$? (x,t, N Q ) and the 
out-of-condensate excitations, e.g., 5rik = (S^l 5ipk) ■ 

Thus, we "dress" the quasistationary parts of the out-of-condensate fields by Eqs. fl54|) 
and (|55| ) and obtain the following set of equations (x — v(x,t)t — > x): 

h 2 - h 2 

ihd t S^ {x,x ± ,t) = \n\5i> (x,t) - — A<5^ (x,t) - — ((v)ipl(x) A<5^o(x,t) + 

+ {(^o + Vo + \fi\ ((v) ) 4> 2 {x) + {2v x + ui) i/>*(x) } 5Mx, x ± , t) + 
+ \v i>Z(x) + 2v x ipo(x) } #o(x, t) + (Jeff ip {x) V5u (x, x ± , t), (56) 
({dt - v(x,t)d x } 2 - cfAj 5u , K (x,x_L,t) - 
- c 2 2k 3 x u o {x) dl5u 0iX (x, x ± , t) - c 2 2k 3 d 2 x u {x) d x 5u 0iX (x, x ± , t) = 

= p' 1 cr eff d x (ip (x) <^ (£,x_L,t) +h.c. ) . (57) 
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Here, we simplified the boson-phonon coupling terms and, roughly, a e g « a e g ~ a . 
Note that the interaction between 5ipQ and 5uo tX is mediated through the condensate 
if) (x, t) ■ d x u Q (x, t) in the linear theory at T — > 0. 

For the outside excitations, the first approximation of Eqs. (BBI) and consists 
of taking uncoupled excitons and phonons. However, they live in a half of the medium 
(x < 0) and form a tail behind the condensate. Then, we have the excitonic excitations 
with 

hu x (k x ) = \fi\ + (k x L ) 2 , \k x \L < 1, 

5ip(x, t) ~ (yV/2 ) _1 exp( i(p k + ik x x — iuu x t ) a^, x < 0. 

This is exactly a free exciton with the energy E g + h 2 ( k — \k x \ ) 2 /2m in the laboratory 
frame. It moves behind the condensate with p x = h(k — \k x \) > 0. For the acoustic 
phonons, we have ik' x > 0): 

hu ph (k' x ) = hu ac (k' x ) -vhk' x , 
5u x (x, xj_, t) ~ const(fc') exp(iip k r + ik' x x — i{u ac (k' x ) — k' x v}t ) b k > + h.c, 



where const(A; / ) oc ( p (V/2) uJac(k') (After x — > x — f t, we have free phonons in 
the tail, and k' x is the same in the laboratory frame.) The 'resonance' equation for the 
excitations written in the comoving frame has the following form: 

\fi\ + \fi\ {k x L ) 2 = h ( Cl T v) k' x , (58) 

and it can be solved easily. Here, we use the following estimate for the 'relevant' phonons: 



(k x L ) w 0.5 



/ mc 1/ 2 
(1 T v/ Cl ) 



For k' x > 0, we have (k' x L ) ~ (1 — 10) y\fi>\ J mcf/2. In this article, we do not go into 

detail of how the coupling between the outside excitons and phonons is formed within 
Eqs. fl56D , (|57[) . Our aim is to show that it is possible to balance the energy of the total 



packet by taking into account the growing out-of-condensate tail ( 51 ) behind the coherent 
state (pi), (O) with leakage (M). 



Returning to Eq. (|53D , we estimate the characteristic energies of excitations as follows: 

{hiv x } ~ h 2 k 2 /2m and {fioJph) ~ const/; const^ ~ 2 — 8. (59) 

(The value of (hup^) comes from the resonance condition at k p ^ x > 0.) We assume 
that the tail consists of excitons and phonons, and the selfconsistency condition between 
ip 2 {x/Lo) oc <3>q and, roughly, d x u o (x/L ) oc $ 2 leads to the coherent emission of the 
outside excitons and phonons, and $ 2 — > $ 2 — S^l(t). Then, we have 

(6iV X)tai i) ^ (SN phMil ). 
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As a result, the energy of the total packet 



E(t) « e Q N Q (t) + 5e N (t) (5v top (t) t/L (t)) 2 -et const iV (t) (5v top (t) t / L (t)f + 

+ (hu x ) const' N (t) (5v top (t) t / L {t)f + (hco ph ) const' N {t) (Sv top (t) t / L (t)) 2 (60) 

is conserved if one takes eo const ~ 5e and const ~ const' ~ 0.1 in Eqs. (|52]) , (|53|) , and 
( pO] ) (the last value is probably overestimated). 

To address the problem of the dynamic stability, one has to estimate how many par- 
ticles and phonons are emitted out from the moving coherent packet. We start from 
the same initial conditions as in Fig. 2. However, the occupations of the tail states 
start growing immediately as t > while the occupancy of the coherent state starts de- 
creasing. Using the conditions E(t) = E Q (t = 0) and N (t) + 5N ta i\(t) = N Q (t = 0), 
we obtain the dynamics presented on Fig. 3 by combination of numerics and analytic 
estimates. For example, after the transport has been observed during the time inter- 
val of At corresponding to 5v top At / L — (1/4 — 1/2), the soliton develops the leak- 
age of 5N ~ (1—3) 10~ 2 N a into the tail excitations, see Fig. 3. However, the value 
of Sv t op{N ) At I L (N o ) starts to decrease because of this leakage too. For instance, 
when the adjusted value of Sv top ( N Q (t) ) At' / L (N o (t)) ~ 1, the leakage is esimated as 
5N ~ (1 — 2) 10 _1 N , see Fig. 3. Therefore, the condensate can be considered as a 
quasistable one during the time interval of ~ At'. The following formulas can be used 
for estimates 

<S>l{t)S ± 2L {t)~N o (t), (61) 

Sn tail {t)S ± L~q{t)N , (62) 

where <5n tail (t) ~ q'§l(t), L = q~L . If q(t) <^ 1 (e.g., q < 10 _1 within the time interval 
At' on Fig. 3) and g > 1 (e.g., q ~ 10 — 50), one can roughly estimate the dimensionless 
q' as ~ lO^ 1 q. The dynamics presented on Fig. 3 is in a qualitative agreement with 
experimental data [[HJ . (It is interesting to mention that the similar problems of dynamic 



stability appear in the theory of the so-called Embedded Solitons in nonlinear optics 
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Figure 3: To model transport properties of the boson-phonon soliton in the case of effec- 
tive dissipation (the "leakage"), we started from the symmetric soliton without a tail as an 
initial condition at t — 0. Dynamics of the boson (exciton) part of the packet is presented 
on Fig. 3 (a) in the form of moving \^ (x — v (x, t)t)\ 2 + 5n(x,t), where 8n(x, t) ~ 5n ta o 
at x < — (2 — 3) L . The coherent phonon part (a moving kink of the displacement field) 
is depicted on Fig. 3 (b); the phonon part of the tail {{d x Su x ) 2 )(t) ^ is not presented 
on this figure. The initial value of the interaction parameter ((v ) $^ ~ Sv top /v is taken to 
be +0.05. Then, the visible changements occur after the packet has traveled the distance 
of (20 — 30) L Q , which corresponds to the effective value of Sv top (t) At/L (t) ~ 1. Note 
that the energy of the total moving packet (coherent part + tail) is conserved at T — > 0. 
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7 Conclusion 



When a quasistationary exciton-phonon soliton moves in a periodic medium, an asym- 
metric form of the ballistic signal can be developed as a dynamic effect due to exciton- 
phonon interaction. This effect has the tendency to be accumulated with time, and it 
can be clearly seen after the packet has traveled the distance of nLo, n ^> 1. Here, the 
width of the packet is ~ (3 — 6) L , and the value of n depends on the strength of relevant 
interaction parameters, such as $o an d ctq- It can be estimated by the following formula 



(see Eq. (@D): 

r2~(|Ho P |/^) _1 ~ 10 3 - 10 2 . 

To a first approximation, we introduced the coherent part of the packet and described 
its dynamics by use of the selfconsistent generalization of v — const — > v(x,t) and 
<p c (x,t) = k x — u t — > (p c (x,t) in such a way: 

v ^ v( d x u (x, *))—«( |*o(:r, t) | 2 ) 

(p c {x,t) ^-2 c {x,t, d x u (x,t)) — >ip c {x,t, \^ (x,t)\ 2 ). 

This substitute was applied within the quasistationary approximation to simplify and solve 
the corresponding dynamic equations on the coherent Bose-core of the total packet. The 
incoherent part of it turned out to be equally important in dynamics of the total packet. 
Indeed, the quasistationary solution for the Bose-core cannot move in a periodic medium 
as a single soliton, but with the leakage into the out-of-condensate excitation states. 
Thus, the question on whether a moving conserving solution (i.e., a soliton conserving at 
least E a , P 0tX , and N Q ) exists in the model with $ 7^ and itq 7^ remains open. In 
this article, we showed that an exciton-phonon "comet" is a better image for the moving 
ballistic packet, and the core of it (a kind of the "nucleus" of such a comet) can be modeled 
by a coherent state, or a condensate. 

The ansatz we used for the soliton-like solution can be applied within the validity of 
the adiabatic hypothesis (i.e., it is not the exact one of the dynamic equations). As a 
result, we could obtain an explicit dependence on time ( cx (Sv top t / L ) 2 ) for the exact 
energy and momentum of the moving soliton if it were the only one component of the 
moving packet. However, the technique of Bose-Einstein condensation allows to proceed 
with the approximate solution obtained in this article. It is known that the interaction 
between Bose-condensate and non-condensed particles leads to an effective dissipation for 



the condensate wave function at T ^ |p0 |. We assume that an effective dissipation 



appears naturally in description of the transport properties of the coherent part, ^ (x,t) 
and d x u Q (x,t), at T — > 0. If the initial state of the packet was taken as a pure coherent 
state with a nonzero momentum, the moving condensate emits excitations, i.e., N Q — > 
N Q — 5N(t) during the observation time. Then, the out-of-condensate cloud of collective 
excitations and, in particular, the tail consisting of the out-of-condensate excitons and 
phonons grow around and behind the moving coherent state, e.g., (5n(x,t) ) 7^ 0, x < 
-3L . 

We argue that the emission effect can make the dynamics of the total moving packet 
conservative. Therefore, the coherent part of such a packet has to be considered as 
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a quasistable core (a "nucleus" of the "comet") of the moving exciton-phonon droplet 
provided 5N ta i\(t) <C N Q (t) during the observation time. In addition, the coherent phase 
<f c (x, t) prescribed to the macroscopic wave function ty Q (x, t) cannot be taken as a regular 
field under such conditions because of dephasing effects Its fluctuations has to be 
taken into account to clarify the coherent properties of the Bose-core. On the other 
hand, the rigorous approach to the dynamics of solitons with emission lies beyond the 
quasistationary approximation we used in this articles. A set of kinetic equations has to 
be applied to treat this problem in detail. 
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